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ABSTRACT 

We report on high-speed eclipse photometry of the dwarf nova V2051 Oph while 
it was in a low brightness state, at B ~ 16.2 mag. In comparison to the average IUE 
spectra, the ultraviolet continuum and emission lines appear reduced by factors of, 
respectively, ~ 4 and ~ 5. Flickering activity is mostly suppressed and the lightcurve 
shows the eclipse of a compact white dwarf at disc centre which contributes ~ 60 per 
cent of the total light at 3900-4300 A. 

We use measurements of contact phases in the eclipse lightcurve to derive the 
binary geometry and to estimate masses and relevant dimensions. We find a mass 
ratio of q = 0.19 ± 0.03 and an inclination of i = 83 ± 2 degrees. The masses of the 
component stars are Mi = 0.78±0.06 Af and M 2 = 0.15±0.03 M . Our photometric 
model predicts K\ = 83 ± 12 km s _1 and K 2 = 435 ± 11 km s _1 . The predicted value 
of Ki is in accordance with the velocity amplitude obtained from the emission lines 
after a correction for asymmetric line emission in the disc is made (Watts et al. 1986). 

The secondary of V2051 Oph is significantly more massive than the secondaries 
of the other ultra-short period dwarf novae. V2051 Oph is probably a relatively young 
system, whose secondary star had not enough time to evolve out of thermal equilib- 
rium. 

Key words: binaries: close - novae, cataclysmic variables - eclipses - stars: individ- 
ual: (V2051 Ophiuchi). 



1 INTRODUCTION 

Cataclysmic Variables (CVs) are close binary systems in 
which mass transfer occurs from a late-type star filling its 
Roche lobe to a companion white dwarf via an accretion 
disc or column (Warner 1996). The subclass of the dwarf 
novae comprises systems in which the white dwarf is non- 
magnetic (B ^ 10 4 G) and which are mostly known for their 
recurrent, large amplitude outbursts (Am = 3 — 5 mag, 
typical duration of 5-10 days). These outbursts reflects a 
change in the structure of their discs - from a cool, optically 



* Based on observations with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-2655, and on ob- 
servations made at Laboratorio Nacional de Astrofi'sica/CNPq, 
Brazil. 



thin, low viscosity state to a hot, optically thick, high vis- 
cosity state - and which are usually parameterized as a large 
change in the mass accretion rate ( M= 10 -11 Mq yr^ 1 i— > 
1CT 9 M Q yr' 1 . See, e.g. Pringle, Verbunt & Wade 1986). In 
quiescence, when the accretion disc is faint, the white dwarf 
and the bright spot (where the stream of transferred matter 
hits the edge of the accretion disc) dominate the optical and 
ultraviolet light from the system, and can be seen as sharp 
changes in the slope of the lightcurve in eclipsing systems. 
From the timings of these features, a precise determination 
of the binary parameters can be obtained. 

V2051 Oph is an eclipsing dwarf nova with an orbital 
period of 90 minutes. High speed photometry by Warner 
(1983) and Cook & Brunt (1983) shows that the orbital 
lightcurve is characterized by deep eclipses (AS ~ 2.5 mag) 
and large amplitude flickering ( ^ 30 per cent), responsible 
for a variety of eclipse morphologies. Significant (> 50 per 
cent) brightness changes occur on timescales shorter than 
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the orbital period but no stable orbital hump (associated to 
anisotropic emission from a bright spot) precedes the eclipse. 
Unlike most dwarf novae, which usually display an orbital 
hump with maximum at phase ~ 0.8, in V2051 Oph the 
humps apparently occur at any orbital phase (Berriman, 
Kenyon & Bailey 1986) . The compact source usually seen at 
disc centre is variable both in brightness and in the dura- 
tion of its immergence and emergence (from ~ 50 to 80 s), 
and was interpreted by Warner & Cropper (1983) to be the 
inner region of the disc and not the white dwarf. Only one 
relatively short outburst of V2051 Oph has been reported 
since its discovery, in which it reached B ~ 13 mag (Warner 
& O'Donoghue 1987). The analysis of optical and infrared 
lightcurves led Berriman et al. (1986) to suggest that the ac- 
cretion disc in V2051 Oph consists largely of optically thin 
material, in accordance with the observed color changes dur- 
ing eclipse (Cook & Brunt 1983). V2051 Oph was proposed 
to be a low field polar system by Warner & O'Donoghue 
(1987) - a possibility first mentioned by Bond & Wagner 
(1977) - based on the interpretation of their eclipse maps 
and on the observation of a 42 s oscillation in the optical 
during outburst, reminiscent of the rapid oscillations seen 
in polars. However, in quiescence the emission lines are dou- 
ble peaked and exhibit the classical rotational disturbance 
effect during eclipse, suggesting the existence of a prograde 
accretion disc (Cook & Brunt 1983; Watts et al. 1986). Watts 
et al. (1986) combined the measured Kj value with velocity- 
dependent eclipse data of the emission lines to derive a dy- 
namical solution for the binary, obtaining Mi = 0.43 Mq, 
M 2 = 0.11 Mq and a radius R 2 = 0.15 Rq. They also fit the 
ultraviolet to optical spectrum with a steady-state optically 
thick disc model to find M= 1.3 x 10 -10 M Q yr' 1 and a 
distance to the system between 90 and 150 pc. 

V2051 Oph belongs to a small group of ultra-short pe- 
riod eclipsing dwarf nova together with Z Cha, OY Car and 
HT Cas. Binary parameters for those other systems have 
been determined with high precision using the method de- 
scribed by Cook & Warner (1984) to derive the mass ratio, 
inclination and white dwarf radius from the white dwarf and 
bright spot eclipse phases. In V2051 Oph, the absence of a 
well defined bright spot and the uncertainty in associating 
the compact central source with the white dwarf (due to 
its variability both in width and brightness) has, until now, 
precluded the use of this method to derive its binary param- 
eters. 

In this paper we present and discuss Hubble Space 
Telescope and ground-based eclipse observations while 
V2051 Oph was in an occasional low brightness state. The 
data are described in section^ In section^] we use contact 
phases of the white dwarf and bright spot to develop a purely 
photometric model for the binary, deriving masses and radii 
of both stars, the binary inclination and the orbital sepa- 
ration. Section^ discusses possible systematic effects affect- 
ing the derived binary parameters, possible causes for the 
observed low brightness, the secondary star and the evolu- 
tionary stage of V2051 Oph. Our results are summarized in 
section |B|. 



2 OBSERVATIONS 



2.1 HST/FOS high-speed spectroscopy 

The Faint Object Spectrograph (FOS) onboard the Hubble 
Space Telescope (HST) was used to obtain time-resolved 
spectroscopy of two consecutive eclipses of the dwarf nova 
V2051 Oph on 1996 January 29. The observations, summa- 
rized in Table jj], consist of two series of exposures in the 
RAPID readout mode at a time resolution of 3.38 s and a 
dead-time between exposures of 0.1s. During data acqui- 
sition, the spectrum was shifted electronically by 1/2 of a 
diode along the 516 diode Digicon array ('1/2 sub-stepping') 
resulting in a 1032-pixel spectrum with a net exposure time 
of 1.64 s pixel -1 . The first run (HI) was carried out using 
a 0.86" circular aperture and the G400H grating (spectral 
resolution AA = 1.5 A pixel -1 ) and covered the egress of 
the first eclipse. The second run (H2) was performed us- 
ing a 0.86" x 0.86" square aperture and the G160L grating 
(AA = 3.5 A pixel -1 ) and covered the following eclipse. The 
spectral range covered by the HST observations is, respec- 
tively for the G160L and the G400H data, 1150 - 2500 A 
and 3200 - 4800 A. 

The observations were reduced using procedures simi- 
lar to the standard STSDAS pipeline, and included flat-field 
and geomagnetically induced motion ('GIMP') corrections, 
background and scattered light subtraction, wavelength and 
absolute flux calibrations. The absolute photometric accu- 
racy of these observations should be about 4 per cent. 

For the G160L observations, the undispersed order-zero 
light is also recorded on the diode array, providing a broad- 
band UV-optical photometric measurement that is simulta- 
neous with each of the ultraviolet spectra. The order-zero 
light has a passband with a FWHM of 1900 A and a pivot 
wavelength of 3400 A. To convert counts to flux densities we 
used the predicted post-COSTAR response of 950 counts s -1 
mjy -1 of Eracleous & Home (1994), which has un uncer- 
tainty of 50 per cent. 

Average out-of-eclipse and mid-eclipse spectra of the 
G160L data are shown in the upper panel of Fig. The 
most prominent lines and bands are labeled. The compari- 
son of the G160L out-of-eclipse spectrum with the average 
IUE spectrum of Watts et al. (1986) reveals that V2051 Oph 
was considerably fainter than usual at the time of the HST 
observations. The flux in the continuum and the emission 
lines have been reduced by factors of, respectively, ~ 4 and 
~ 5. The lower panel of Fig[l] shows average out-of-eclipse 
and mid-eclipse spectra of the G400H data. The G160L 
order-zero passband overlaps the spectral coverage of the 
G400H data. Filled squares show the fluxes of the G160L 
order zero light computed at the same orbital phase ranges 
as the G400H average spectra. This comparison indicates 
that V2051 Oph remained at roughly the same brightness 
level during the HST observations. 

Lightcurves at selected passbands for both HST runs 
are shown in Fig.^. These will be discussed in detail in sec- 
tion 2.4. 



2.2 high-speed photometry 

High-speed photometry of V2051 Oph was obtained with 
the 1.6 m telescope of Laboratorio Nacional de Astrofi'sica 
(LNA/CNPq), in southern Brazil, on 1996 June 15-17. Four 
eclipses were collected, one in the B-band and the remain- 
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Table 1. Journal of the observations. 



date 


run 


Start 


At 


No. of 


Spectral 


Phase range 


Instrument 


(1996) 




(UT) 


(s) 


exposures 


range 


(cycles) 




Jan 29 


HI 


19:03 


3.38 


113 


3226-4781 A 


+0.01, +0.08 


HST/FOS/BL 


Jan 29 


H2 


20:24 


3.38 


693 


1150-2507 A 


-0.09, +0.34 


HST/FOS/BL 


Jun 15 


LI 


22:21 


5.0 


547 


B-band 


-0.24, +0.27 


FOTRAP 


Jun 16 


L2 


22:24 


5.0 


398 


white light 


-0.18, +0.18 


FOTRAP 


Jun 16 


L3 


23:45 


5.0 


516 


white light 


-0.29, +0.19 


FOTRAP 


Jun 17 


L4 


22:19 


5.0 


432 


white light 


-0.23, +0.17 


FOTRAP 



ing in white light (W). These observations are detailed in 
Table [j]. 

The observations were performed with the one-channel 
FOTRAP photometer in its single-filter mode at a time reso- 
lution of 5 s, using a 11.3" diaphragm. A Hamamatsu R943- 
02 PMT was employed. Sky measurements were taken at 
intervals of 15-20 min at a position 20" to the north of the 
variable, except during eclipse. The sky measurements were 
fitted by a low-order polynomial and then subtracted from 
the raw data. A close comparison star was also regularly 
observed to check for the presence of clouds and sky trans- 
parency variations. From these observations, we confirmed 
that all 4 runs were performed under good sky conditions. 
Coefficients of extinction and of transformation to the UB- 
VRI standard system were derived from observations of E- 
regions stars of Graham (1982) and blue spectrophotometric 
standards of Stone & Baldwin (1983). The reader is referred 
to Jablonski et al. (1994) for a detailed study of the reli- 
ability of the transformations from FOTRAP's natural to 
the standard UBVRI system, in particular for objects of pe- 
culiar spectra such as CVs. We used the relation of Lamia 
(1981), B = 2.51og/ s (mJy) - 16.57, to transform B magni- 
tudes to flux density units. The absolute calibration of the 
B-band observations is accurate to better than 4 per cent. 

The ground-based lightcurves of V2051 Oph are shown 
in Fig.W and will be discussed in section 2.4. 



2.3 New eclipse timings 

Mid-eclipse times of V2051 Oph were measured by employ- 
ing the bissected chord method and by computing the mid- 
point of the phases of maximum and minimum derivative in 
the lightcurves (section^). These methods yield values which 
are consistent with each other under the uncertainties, and 
therefore we adopted the mean as the mid-eclipse time. The 
new heliocentric (HJD) as well as baricentric (BJED) tim- 
ings of V2051 Oph are listed in Table^J with corresponding 
cycle number and uncertainties (quoted in parenthesis). Also 
shown are the (O-C) values with respect to the revised lin- 
ear ephemeris of Baptista, Tripplet & Bond (1998), 

T mid = BJED 2 443 245.977 850 + 0.062 427 8595 x E (1) 

The measured (O-C) values are in accordance with the cycli- 
cal variations in the orbital period of V2051 Oph observed 
by Baptista, Tripplet & Bond (1998). 



2.4 Eclipse lightcurves 

We adopted the following convention regarding the phases: 
conjunction occurs at phase zero, the phases are negative 
before conjunction and positive afterwards. The lightcurves 



Table 2. 


New eclipse timings. 




Cycle 


HJD 


BJED 


(o-c)R 




(2,450,000+) 


(2,450,000+) 


(cycles) 


109 988 


112.29277 


112.29273(4) 


-0.0084 


109 989 


112.35517 


112.35513(4) 


-0.0089 


112 201 


250.44558 


250.44555(2) 


-0.0090 


112 217 


251.44443 


251.44440(1) 


-0.0089 


112 218 


251.50688 


251.50685(1) 


-0.0086 


112 233 


252.44331 


252.44328(2) 


-0.0084 



a Observed minus calculated times with respect to the 
ephemeris of eq.JI]. 



were phased according to the linear ephemeris of eq. |jj A 
phase correction of +0.0087 cycles was further applied to the 
data to make the centre of the white dwarf eclipse coincident 
with phase zero (see section^). 

V2051 Oph was particularly faint at the time of both 
the HST and the ground-based observations. From the B- 
band lightcurve, we estimate a magnitude of B = 16.3 mag, 
while the G400H data yield B ~ 16.2 mag. This is consider- 
ably fainter than the B = 14.7 — 15.5 mag range previously 
reported in the literature (Warner & Cropper 1983; Cook & 
Brunt 1983; Wenzel 1984; Berriman et al. 1986; and Watts 
et al. 1986). We will thereafter refer to this as the 'low' or 
'minimum' brightness state of V2051 Oph to differentiate it 
from the normal, quiescent state. Thus, our data yield a rare 
(and lucky) opportunity to cleanly see the white dwarf and 
bright spot in V2051 Oph, with reduced contributions from 
the accretion disc and boundary layer to the total light. 

Figure^ shows HST lightcurves at selected passbands, 
wavelength range increasing upwards. The eclipse shape is 
reminiscent of those of the dwarf novae OY Car and Z Cha 
in quiescence (e.g., Wood et al. 1986, 1989), with two com- 
pact sources being eclipsed in sequence. Sharp steps corre- 
sponding to the egress of the white dwarf and bright spot are 
clearly seen in the G400H lightcurves. A substantial fraction 
of the emitted light comes from the white dwarf, which con- 
tributes ~ 60 per cent of the total light at 3900-4300 A. The 
G160L order zero lightcurve shows the sequence of ingress 
of two compact sources. The phases of ingress of the first 
eclipsed source are coincident with those of the white dwarf 
as estimated from the ground-based data (Fig^). Hence, we 
can unambiguously associate the first compact source with 
the white dwarf and, therefore, the source which is eclipsed 
later is the bright spot. Egress features of both sources can 
also be easily identified in the order zero lightcurve. The flux 
of the bright spot is larger at ingress than at egress. This 
is consistent with the overall decrease in brightness along 
the run and is usually understood in terms of anisotropic 
emission due to the foreshortening of a spot 'painted' at the 
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edge of an optically thick accretion disc. White dwarf and 
bright spot features are less well defined in the other, lower 
signal-to-noise ratio G160L lightcurves. Flickering activity - 
usually higher in the ultraviolet than at optical wavelengths 
in CVs - is considerably less impressive than observed in 
quiescence (e.g, Warner & O'Donoghue 1987), probably re- 
flecting the lower brightness state of the system at the epoch 
of our observations. 

The HST dataset gives a clear and simultaneous mea- 
surement of the white dwarf and bright spot eclipse phases, 
which will allow us to estimate the binary mass ratio in sec- 
tion |. 

Figure^ display the ground-based observations. These 
lightcurves clearly show the total eclipse of a compact source 
at disc centre, the white dwarf, and reveal that the white 
dwarf contributes 80-90 per cent of the light at that epoch, 
with very little evidence of an accretion disc or bright spot. 
Flickering activity is mostly suppressed, which is possibly 
an indication that mass transfer was substantially reduced 
at that epoch. The duration of the ingress/egress feature 
(~ 30 s or 0.006 cycles, see section^) is shorter than that in 
quiescence (Warner & Cropper 1983; Berriman et al. 1986; 
Watts et al. 1986) by a factor of two, possibly an indication 
that, at those epochs, one was seeing a bright (and larger) 
boundary layer around the white dwarf instead of the white 
dwarf itself. This possibility was previously mentioned by 
Warner & Cropper (1983). 

The ground-based dataset yields an unambiguous mea- 
surement of the total width of the eclipse of the white dwarf 
as well as of the duration of its ingress/egress feature and 
will be crucial for our derivation of the binary parameters 
in section 



3 THE BINARY PARAMETERS 

In this section we use measurements of the contact phases of 
the white dwarf and of the bright spot to develop a model for 
the binary in V2051 Oph based solely on quantities obtained 
from the photometry. 

3.1 Measuring contact phases 

The ingress/egress phases of the occultation of the white 
dwarf (hereafter WD) and of the bright spot (BS) by the 
secondary star provide information about the geometry of 
the binary system and the relative sizes of these components 
(Wood et al. 1986, 1989). The contact phases can be iden- 
tified as rapid changes in the slope of the lightcurve, visible 
in the lightcurves shown in Figs. ^ and ^j. We follow the 
notation of Baptista et al. (1989), taking w i,0w2 as those 
phases during which WD disappears behind the secondary 
star and d> the phases corresponding to the reap- 

pearance of WD from eclipse. We similarly define a set of 
phases 0t>i-4 for BS. 

The contact phases were measured with the aid of a 
cursor on a graphic display of a median filtered version of 
the lightcurve. The estimated error for this procedure is 
0.0006 cycles. Mid-ingress (egress) phases for WD and BS 
(c/>wi, <t>uDe and cf>bi,4>be) were computed as the mid point of 
the ingress (egress) feature and also as the phase at which 
half of the compact source light is eclipsed. 



We also employed the technique described by Wood, 
Irwin & Pringle (1985) to measure the contact points and 
to estimate the phases of mid-ingress (egress), here defined 
as the points of minimum (maximum) derivative in the 
lightcurve. Figure^] illustrates the procedure for the G160L 
order zero lightcurve. The original lightcurve [Fig.^a)] is 
smoothed with a median filter of width 0.0019 cycles (3 data 
points) [Fig.^(b)] and its numerical derivative is calculated. 
The amount of filtering used at this point is a compromise 
between the aim to suppress noise in the lightcurve as much 
as possible and the need to preserve any real structures in 
the eclipse shape. The ingress/egress of WD and BS can be 
seen as those intervals for which the derivative curve is sig- 
nificantly different from zero. The width at half-peak inten- 
sity of these features yields a preliminary estimate of their 
duration. A spline function is fitted to the remaining regions 
in the derivative curve to remove the contribution from the 
extended and slowly varying eclipse of the disc. The spline- 
subtracted derivative curve [Fig.W(c)] is then analyzed by an 
algorithm which identifies the points of extrema (the mid- 
ingress / egress phases) and the points where the derivative 
starts to depart from the zero-constant level (the contact 
points). A median filter of equal width to that used in the 
first part of the reduction is applied to the derivative curve 
to aid in the detection. A similar sequence is used to measure 
the contact points and mid phases of BS. 

Table g collects the measured contact phases, mid- 
ingress and mid-egress phases of WD. Quoted values are the 
average of the determinations from all methods described 
above. Uncertain measurements are marked by a colon. 
Figja depicts the median contact phases of WD listed in Ta- 
bleHin an enlarged view of the combined W+B lightcurve 
and its derivative curve. 

Integrating the flux in the derivative curve Fig.^(c) be- 
tween the first and second (third and fourth) contact phases 
we obtain estimates of the WD flux at ingress (egress). The 
lightcurve of the WD is then reconstructed by assuming that 
the flux is zero between ingress and egress and that it is con- 
stant outside eclipse. The reconstructed WD lightcurve can 
be seen in Fig. ^(d) , and the lightcurve after removal of the 
WD component is shown in Fig.^e). 

For the HST order zero lightcurve (fig. ^|), the estimated 
WD flux at ingress is 0.40 ±0.05 mjy, larger than the flux of 
0.30 ±0.05 mjy obtained from the egress feature. This could 
be due to a phase dependent absorption of the WD light, 
the absorption being more pronounced at egress phases. The 
analysis of Catalan et al. (1998) suggests that the WD in 
V2051 Oph is veiled by a large number of blended Fell lines 
(the so called Fell curtain). This effect is reminiscent of that 
observed previously in OY Car (Home et al. 1994) and has 
been attributed to absorption by circumstellar material, pos- 
sibly in the outer disc. The observed difference in WD flux 
at ingress and at egress could be associated with this effect if 
the absorption by the Fe II curtain is more pronounced along 
the line of sight at egress phases. This seems to be the case 
in the novalike UX UMa, whose ultraviolet eclipse mapping 
reveals that the absorption by the Fell curtain in the disc 
side closest to the secondary star is more pronounced in the 
leading quadrant than in the trailing quadrant (which con- 
tains the bright spot and gas stream) and that this effect is 
probably responsible for a substantial reduction in flux level 
from before to after eclipse (Baptista et al. 1998b). 
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Table 3. Contact phases of the white dwarf. 



run 


spectral range 


4>wi 


4 l w2 


<t>w3 


<t>w4 




4* we 


HI 


3236-3600 A 


- 


- 


+0.0312 


±0.0369 


- 


±0.0339 




3900-4300 A 


- 


- 


±0.0303 


±0.0365 


- 


±0.0335 




4400-4781 A 


- 


- 


±0.0309 


±0.0363 


- 


±0.0337 




3236-4781 A 


- 


- 


±0.0306 


±0.0361 


- 


±0.0337 


H2 


1400-1800 A 


-0.0348 


-0.0317: 


±0.0285 


±0.0347 


-0.0332 


±0.0331 




1802-2303 A 


-0.0361 


-0.0298 


±0.0279 


±0.0361 


-0.0330 


±0.0325 




1152-2506 A 


-0.0359 


-0.0303 


±0.0280 


±0.0341 


-0.0328 


±0.0334 




order zero 


-0.0366 


-0.0309 


±0.0305 


±0.0355 


-0.0334 


±0.0330 


LI 


B 


-0.0374 


-0.0309 


±0.0307 


±0.0367 


-0.0332 


±0.0332 


L2 


W 


-0.0359 


-0.0285 


±0.0308 


±0.0369 


-0.0322 


±0.0336 


L3 


W 


-0.0356 


-0.0300 


±0.0302 


±0.0367 


-0.0333 


±0.0334 


L4 


w 


-0.0354 


-0.0303 


±0.0304 


±0.0369 


-0.0328 


±0.0332 


combined W+B 




-0.0355 


-0.0298 


±0.0306 


±0.0372 


-0.0328 


±0.0334 


median 




-0.0359 


-0.0301 


±0.0304 


±0.0364 


-0.0330 


±0.0334 


error 




±0.0004 


±0.0005 


±0.0003 


±0.0003 


±0.0002 


±0.0002 



Table 4. Contact phases of the bright spot, 

run spectral range <j>i,i <f>b2 

HI 3236-3600 
3900-4300 
4400-4781 
3236-4781 

H2 order zero -0.0295 -0.0226 



4>b3 4>bi <Pbi 4>be 

±0.0811 ±0.0855 ±0.0831 

±0.0804 +0.0848 ±0.0822 

±0.0811 ±0.0836 ±0.0825 

±0.0809 ±0.0831 ±0.0821 



±0.0841 +0.0892 -0.0260 +0.0866 



Measured BS contact phases, mid-ingress and mid- 
egress phases are listed in Table ^. As before, quoted val- 
ues are the average of the determinations from all meth- 
ods described above. For run H2, only values for the order 
zero lightcurve are listed since the lightcurves at the other 
passbands (Fig.|^) do not yield reliable measurements of BS 
phases. The ingress phases fai, 4>b2 and (f>u were measured 
from the WD subtrated lightcurve (Fig.^(e)), which pro- 
vides an unblended, clean view of the BS ingress feature. 
The short run HI gives only the egress phase of the BS, 
which occurs perceptibly earlier than at run H2. Warner 
& Cropper (1983) and Warner & O'Donoghue (1987) have 
shown that the eclipse shape of V2051 Oph (and therefore its 
disc structure) is very variable and that significant changes 
occur on the timescale of one orbit. One expects that sim- 
ilar changes in disc radius and BS position do occur in the 
same timescale. We have thus been cautious not to combine 
egress phases from different cycles and, theref ore, for our 
procedure to determine the mass ratio q (section 3.2) we will 
use only the pair of mid-ingress/egress phases of BS mea- 
sured from run H2. The measurements of mid-ingress/egress 
phases from the three different methods have dispersions of 
0.0005 cycles. However, since only one light curve was used 
to estimate these values, we adopted more conservative er- 
ror bars of 0.001 cycles for the mid-ingress phase and 0.002 
cycles for the mid-egress phase. 

The duration of the WD eclipse is defined as 



'we ywi i 



(2) 



and the centre of the WD eclipse (the inferior conjunction 
of the binary) is written as 

4>o — 1/2 (4> we + 4> W i) ■ (3) 
These quantities are collected in Table |^. The median of 



the measurements from all lightcurves yields A(f> = 0.0662 + 
0.0002 cycles, where the quoted error is the median of the 
absolute deviations with respect to the median. Similarly, 
we have <j>o — +0.0000 ± 0.0002 cycles, which indicates that 
the centre of the WD eclipse corresponds to phase zero. 

The difference between the first and second (third and 
fourth) WD contact phases yield the phase width of the 
WD ingress (egress), A w i (A wo ). These quantities are also 
listed in Table ^| A median from all values of A w i and A wc 
yields A wt j = 0.0060 ± 0.0005 cycles, where again the error 
is the median of the absolute deviations with respect to the 
median. 



3.2 Mass ratio, inclination and disc radius 

Making the usual assumption that the secondary star fills 
its Roche lobe and given the duration of the eclipse of the 
central parts of the disc, A</>, there is a unique relation be- 
tween the mass ratio q = M2 jM\ and the binary inclination 
i (Bailey 1979; Home 1985). From Table|, the width of the 
eclipse in V2051 Oph is A<f> = 0.0662. This gives the con- 
straint q > 0.116, with q = 0.116 if i = 90°. 

When combined with the measured contact phases of 
the WD and BS, this relation gives a unique solution for q, 
i, and Rbs/RLi, where Rbs is the distance from disc centre 
to the BS (usually taken to be the disc radius) and Rli is 
the distance from disc centre to the inner lagrangian point 
LI. The following assumptions are made: (i) the infalling gas 
stream leaving the LI point describes a balistic trajectory 
in the primary lobe (Lubow & Shu 1975), and (ii) the bright 
spot is located at the position at which the stream hits the 
edge of the disc. The theoretical trajectories described by a 
test particle in the primary lobe are solely determined by 
the mass ratio of the binary. For increasing q values the 
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Table 5. White dwarf eclipse parameters. 



run 


spectral range 


A0 


00 


A wi 




HI 


3236-3600 








0.0057 




3900-4300 








0.0062 




4400-4781 








0.0054 




3236-4781 








0.0056 


H2 


1400-1800 


0.0663 


-0.0001 


0.0032: 


0.0062 




1802-2303 


0.0655 


-0.0003 


0.0063 


0.0082 




901-2506 


0.0662 


+0.0003 


0.0056 


0.0061 




order zero 


0.0664 


-0.0002 


0.0057 


0.0050 


LI 


B 


0.0664 


+0.0000 


0.0065 


0.0060 


L2 


W 


0.0658 


+0.0007 


0.0074 


0.0061 


L3 


w 


0.0667 


+0.0000 


0.0056 


0.0065 


L4 


w 


0.0660 


+0.0002 


0.0051 


0.0065 


combined W+B 




0.0662 


+0.0003 


0.0057 


0.0066 


median 




0.0662 


+0.0000 


0.0057 


0.0061 


error 




±0.0002 


±0.0002 


±0.0006 


±0.0007 



trajectories are closer to the line joining both stars, as the 
angular momentum of the infalling gas stream after leaving 
the LI point is progressively smaller. The correct mass ratio, 
and hence inclination, are those for which the calculated 
stream trajectory passes through the observed position of 
the bright spot. This method was first developed by Smak 
(1971) and has since then been applied to various CVs (e.g., 
Fabian et al. 1978; Cook & Warner 1984; Wood et al. 1986, 
1989; Baptista, Steiner & Cieslinski 1994). 

Figure^ illustrates the determination of q, i and 
Rbs/Rll from the eclipse phases. Fig.^(a) shows a diagram 
of ingress versus egress phases for the measurements of the 
WD and BS in Tables| and §. Measurements of the WD 
ingress and egress are shown as the cluster of small crosses 
around phases (—0.033, +0.033) in the lower portion of the 
diagram. The small box in the upper left corner marks the 1- 
a range of the mid-ingress and mid-egress phases of the BS. 
Theoretical gas stream trajectories corresponding to a set of 
pairs (i, q) are also shown. The trajectories were computed 
by solving the equations of motion in a coordinate system 
synchronously rotating with the binary, using a 4th order 
Runge-Kutta algorithm (Press et al. 1986) and conserving 
the Jacobi integral constant to one part in 10 6 . 

The theoretical trajectory that passes through the po- 
sition of the BS has q = 0.19 ± 0.03 and i = 83?3 ± 1?4, 
where the uncertainties are taken from the standard devia- 
tion of the points about the trajectory of best fit. Fig.|^(b) 
shows the geometry of the binary system for q — 0.19. For 
this mass ratio, the relative size of the primary Roche lobe 
is Rli/cl — 0.66, where a is the orbital separation. 

The squashed circle in Fig.^|(a) represents the accre- 
tion disc whose edge passes through the position of BS for 
the adopted mass ratio. This corresponds to a disc radius 
of Rbs/Rhi = (0.56 ± 0.02). A circle with this radius is de- 
picted in Fig.^](b). The calculated radius is halfway between 
the radius expected for zero- viscosity discs, Rd/Rhi =0.27 
(Flannery 1975), and that expected for pressureless discs, 
Rd/RLi=0.76 (Paczyriski 1977). 



3.3 Masses and radii of the component stars 

The duration of the ingress / egress of WD yields an estimate 
of the relative size of the white dwarf. This information can 



be used together with the calculated mass ratio to constraint 
the masses and radii of the component stars. We use the 
approximate relations of Ritter (1980), 



27rA wd ~ 2R 1 /[z(q)a sin 9} , 



and 



Rz 



(4) 



(5) 



where i?i and R2 are the radii of the WD and of the sec- 
ondary star, and z(q) is the relative distance from disc cen- 
tre to the point tangent to the surface of the secondary that 
marks the beginning/end of the eclipse of WD (Baptista et 
al. 1989). The value of z(q) is a by-p roduct of the inclination- 
determination routine (section 3.2) and is usually close to 
unity. In our case, for q — 0.19 we have z = 0.96. The er- 
ror introduced by the use of these expressions is of order of 
(R1/R2) 2 and for the parameters listed in Tabled is smaller 
than 0.5 per cent. 

Combining Eq. (Q) with Kepler's third law we can write 
an expression relating the mass and the radius of the white 
dwarf for V2051 Oph, 



Ri/R@ = 2.349 A wd (Mi/M ch ) 1/3 /(g) , 



(6) 



: 0.0112 



(7) 



where M c h = 1.44 Mq is the Chandrasekhar mass limit and 
f(q) = z sin 8 (1 + q) 1//a . Another relation between Ri and 
Mi can be obtained from Nauenberg's (1972) analytical ap- 
proximation to the mass-radius relation for cool degenerate 
white dwarfs of Hamada & Salpeter (1961), 

iM~ 2/3 „ ( iM 2/31 1/2 

AW W/cJ 

Combining Eqs. (^) and (|?j) one can solve for Mi (q, A w d) 
and the remaining system parameters. 

The primary-secondary mass diagram for V2051 Oph 
can be seen in Fig. [?]. A Monte Carlo propagation code was 
used to estimate the errors in the calculated parameters. 
The values of the input parameters q and A wt j are inde- 
pendently varied according to Gaussian distributions with 
standard deviation equal to the corresponding uncertain- 
ties. The results, together with their 1-cr errors, are listed in 
Table []. The cloud of points in Fig.|^ was obtained from a 
set of 10 4 trials using this code. The highest concentration 
of points indicates the region of most probable solutions. 
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Table 6. The calculated parameters of V2051 Oph. 



Mi/M Q 
Ri/Rq 
M 2 /M Q 
R2/R0 



0.19 ±0.03 
83.3° ± 1.4° 
0.78 ± 0.06 
0.0103 ± 0.0007 
0.15 ±0.03 
0.16 ±0.01 



a/R Q 
Rli/o 

Rbs/RLl 

Ki (km/s) 
K 2 (km/s) 



0.64 ±0.02 
0.66 ±0.01 
0.56 ±0.02 
19° ± 1° 
83 ± 12 
436 ± 11 



4 DISCUSSION 

4.1 The reliability of the photometric model 

Our photometric model is based on the inferred A w< j and q 
values. Here we discuss systematic errors which may affect 
these quantities and their influence on the derived param- 
eters of V2051 Oph, and compare the derived parameters 
with those obtained by other authors. 

Errors in q are related to errors in measuring mid- 
ingress/egress phases, in particular the BS phases. WD 
eclipse phases were precisely measured in the ground-based 
lightcurves, where the BS and the accretion disc give neg- 
ligible contribution to the total light. These phases helped 
us to identify the WD eclipse features in run H2 without 
ambiguities. The lightcurve after the subtraction of the WD 
contribution (fig.^(d)) clearly shows the asymmetric eclipse 
of the compact BS and yielded accurate measurements of 
its mid-ingress/egress phases. The measured BS phases map 
into a small region in the ingress/egress phases diagram from 
which a precise value of q was obtained. Therefore, the mass 
ratio seems well constrained. 

The remaining binary parameters were derived under 
the assumption that the duration of the ingress/egress of the 
central object, A w( j, gives a direct measure of the radius of 
the white dwarf. This supposition may not always be correct. 

It may be possible that the observed compact object of 
Fig. P is a hot and opaque boundary layer involving the white 
dwarf, in which case the true white dwarf would be smaller 
than assumed (Wood & Crawford 1986). In this case, the 
primary mass found is a lower limit. However, notice that 
increasing the white dwarf mass along the line q — 0.19 
in Fig[j] results in a secondary star which is overmassive in 
comparison to a main sequence star of same radius. This 
would be in marked contrast with the other short-period 
eclipsing CVs, where the results seem to point in the oppo- 
site direction, i.e., secondary stars which are undermassive 
in comparison to main sequence stars of same radius (Wood 
et al. 1986, 1989; Home, Wood & Stiening 1991). Other ef- 
fects that may result in an overestimated white dwarf radius 
include the questionable application of a mass-radius rela- 
tion for cool white dwarfs to a relatively hot object (Home, 
Wood & Stiening 1991), and a possible spherical distortion 
for a white dwarf rotating close to breakup velocity (Wood 
& Home 1990). However, the corrections in the derived pa- 
rameters implied by these effects ( ^ 5 per cent) are smaller 
than the uncertainty in the value of A m d and were not taken 
into consideration. 

Alternatively, could the white dwarf be larger than in- 
ferred from our value of A m d! We now address two possibil- 
ities in this regard. 

Previous eclipse photometry of V2051 Oph indicate a 
duration of the ingress/egress of the compact central source 
of A wd ~ 0.009 - 0.012 (Warner & Cropper 1983; Cook & 



Brunt 1983; Watts et al. 1986), substantially larger than 
measured from our observations. However, the variability 
and the clear correlation between brightness and duration of 
ingress/egress of this source led Warner & Cropper (1983) 
to suggest that it might be a bright (and variable) boundary 
layer around the white dwarf and not the white dwarf itself. 
This scenario is supported by our observations at minimum 
light, when mass accretion has probably been substantially 
reduced and the boundary layer has accordingly dimmed, 
leaving a bare, fainter, and much smaller white dwarf as the 
only light source at disc centre. 

Warner & O'Donoghue (1987) suggested that V2051 
Oph could be a low field polar. Following their suggestion, 
we explored the possibility that the observed compact source 
at disc centre is in fact a bright spot at the surface of the 
white dwarf, corresponding to the pole of an accretion col- 
umn. If the white dwarf rotates synchronously with the bi- 
nary, one would expect the lightcurves to show a strong or- 
bital modulation due to the self occultation of the accretion 
column. For example, at phase <f> = 0.25 about half of the 
accretion column would have disappeared behind the white 
dwarf resulting in a reduction of XL 40 per cent in flux of 
the lightcurve. On the other hand, if the white dwarf is non 
synchronously rotating, one would expect to see two conspic- 
uous effects in the lightcurves: (1) large amplitude (~ 40 per 
cent) pulsations with the rotation period of the white dwarf, 
and (2) measurable wanderings of the ingress/egress phases 
and eclipse width of the compact source at disc centre as 
the accretion column is seen at different aspects. Our obser- 
vations clearly exclude both possibilities: there is no sign of 
strong orbital modulation, no sign of pulsations, and no de- 
tected changes in eclipse phases. Therefore, we are left with 
the conclusion that the compact source at disc centre of our 
observations is indeed the white dwarf. 

It is worth mentioning that our observations do not 
discard the possibility that V2051 Oph possess a magnetic 
white dwarf. If the observed minimum brightness state cor- 
responds to a phase of reduced mass transfer and negligible 
or no mass accretion, then no emission would be produced 
in the accretion column, which could therefore remain un- 
detected. 

We now turn to the comparison of our derived param- 
eters with those obtained by Watts et al. (1986). Their dy- 
namical model of V2051 Oph is based on spectroscopy of 
the accretion disc's emission lines. Radial velocity for the 
H/3 and H7 lines yields a K\ = 111 km s _1 and a phase 
lag of +0.13 cycles with respect to the photometric phase. 
Phase lags of 0.1 — 0.2 cycles between photometric and spec- 
troscopic conjunction are seen in many of the CVs for which 
emission lines radial velocity curves were measured. Asym- 
metries in the line intensity distribution (which displace the 
light centroid from disc centre) or the existence of a non- 
circular component to the velocity of the emitting gas (as 
it may happen in the presence of a disc wind) lead to sys- 
tematic errors in the amplitude and phase of the measured 
K\. They show that the observed phase lag of V2051 Oph 
can be explained if the accretion disc is asymmetric in the 
emission lines, with the trailing lune being ~ 60 per cent 
brighter than the leading lune. The measured K\ will be 
overestimated by ~ 22 per cent in this case. The correction 
of their K\ value by this effect eliminates the phase lag and 
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brings their measurement to K\ =91 kms -1 , into good 
accordance to our prediction in Table^. 

Watts et al. (1986) used the double peaked Balmer line 
profiles and their behaviour during eclipse to constrain the 
mass ratio of the binary. Their model predicts a larger mass 
ratio (q = 0.26 ± 0.04) and a much smaller white dwarf 
mass (Afi = 0.43 ± 0.05 M Q ) than ours. We note that their 
fit to the line behaviour during eclipse seems of rather poor 
quality. The larger spectroscopic mass ratio predicts that 
the bright spot eclipse occurs earlier than observed, with a 
phase difference at the 2-a level. The spectroscopic model 
also predicts a duration of the ingress and egress of the white 
dwarf of — 0.011, which clearly conflicts with the ob- 
served duration of A wd = 0.0060 ± 0.0005. 

It would be useful to have another observational con- 
straint on the binary parameters. The best way of checking 
the reliability of the above assumptions in calculating binary 
parameters is to measure the radial velocity of the secondary 
star, K2- 

4.2 The low brightness state: reduced mass 
transfer rate? 

At the time of our observations V2051 Oph was in a much 
lower brightness state than at quiescence, when it has V ~ 
15 mag. Since it was observed to be in this lower brightness 
state both in 1996 January and in 1996 June, it may be 
possible that this is a relatively prolonged minimum state 
of duration a few months. Recovery of photometric observa- 
tions of V2051 Oph during this period would be important 
to clarify this point. 

The system has been observed to go into outburst once, 
where it increased in brightness to reach V ~ 13.2 (Warner 
& O'Donoghue 1987). If the outburst is caused by the disc 
instability mechanism, than the disc is already in a state 
of low viscosity, low mass accretion rate when the system 
is at V ^ 15 mag. Therefore, it is hard to explain the ob- 
served low brightness state in terms of disc instability. The 
simplest way to understand the reduced brightness level of 
V2051 Oph is to consider it the consequence of a significant 
reduction in the mass transfer rate. 



4.3 The secondary star mass and its evolutionary 
consequences 

Our photometric model indicates that V2051 Oph has a 
larger mass ratio and a more massive white dwarf than the 
other eclipsing dwarf novae of similar period. Probably more 
relevant, it has a secondary star which is significantly more 
massive than those of Z Cha, OY Car and HT Cas. Since 
the inferred size of the secondary in these systems is simi- 
lar, R2/RQ — 0.15, this means that the secondary star of 
V2051 oph obeys a mass-radius relationship which predicts 
significantly denser stars for a given mass. 

Theoretical evolutionary models of CVs predict that 
the long term, continuous mass loss brings the secondary 
stars out of thermal equilibrium and make these stars larger 
than a main sequence star of same mass (Rappaport, Joss 
& Webbink 1982). Hence, the fact that the secondary stars 
in Z Cha, OY Car and HT Cas are oversized in compari- 
son to main sequence stars of same mass suggests that these 



systems are relatively old cataclysmic variables, for which 
mass transfer (and loss from the secondary) has been occur- 
ring for some 10 9 years. On the other hand, the fact that 
the secondary star of V2051 Oph obeys a typical main se- 
quence mass-radius relationship suggests that V2051 Oph is 
a relatively young system, for which mass loss has not oc- 
curred for long enough time to bring the secondary star out 
of thermal equilibrium. Further support to this suggestion 
comes from the following argument. If we assume that CVs 
evolve by loosing mass in recurrent nova outbursts (at a rate 
of ~ 10~ 5 Mq at every 10 4 - 10 s years), than V2051 Oph 
will need about 10 9 years to reduce its total mass of about 
0.9 Mq to the typical total mass of ~ 0.7 Mq of the other, 
older dwarf novae of similar orbital period. 



5 CONCLUSIONS 

The analysis of HST/FOS and ground-based eclipse obser- 
vations of the dwarf novae V2051 Oph yielded the following 
results: 

(i) V2051 Oph was observed in a low brightness state, at 
B ~ 16.2 mag. In comparison to the average IUE spectra, 
the ultraviolet continuum and emission lines appear reduced 
by factors of, respectively, ~ 4 and ~ 5. Flickering activity 
is mostly suppressed and the lightcurve shows the eclipse 
of a compact white dwarf at disc centre which contributes 
~ 60 per cent of the total light at 3900-4300 A. In 1996 
January, the bright spot is clearly seen in the lightcurve, 
allowing a precise measurement of its eclipse phases. In 1996 
June only the white dwarf is seen in the lightcurve and its 
radius can be estimated from accurate measurement of the 
width of its ingress/egress. It is suggested that this state is 
the consequence of a substantial reduction in mass transfer 
rate. 

(ii) We developed a purely photometric model for the bi- 
nary from the eclipse phases of the white dwarf and bright 
spot. For the gas stream trajectory to pass through the po- 
sition of the bright spot a mass ratio q — 0.19 ± 0.03 is 
required. The corresponding inclination is i — 83?3 ± 1?4. 

(iii) The duration of the ingress/egress of the white dwarf 
combined with the Hamada-Salpeter mass-radius relation- 
ship and the obtained mass ratio give Mi = 0.78 ± 0.06 Mq 
and M2 = 0.15 ± 0.03 Mq. Calculated parameters are sum- 
marized in Table [| 

(iv) The secondary of V2051 Oph is significantly more 
massive than the secondaries of the other ultra-short period 
dwarf novae and seems to obey a ZAMS mass-radius rela- 
tionship. V2051 Oph is probably a relatively young system, 
whose secondary star had not enough time to evolve out of 
thermal equilibrium. 

(v) Our photometric model predicts K\ — 83 ± 12 km s _ , 
which is in accordance with the velocity amplitude obtained 
from the emission lines after a correction for asymmetric line 
emission in the disc is made (Watts et al. 1986). 
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Figure 1. Average out-of-eclipse (black) and mid-eclipse (light gray) HST spectra of V2051 Oph. Major emission and absorption features 
are labelled. Phase ranges of each spectrum are as follows: +0.1 to +0.3 (G160L, out of eclipse), —0.02 to +0.03 (G160L, mid-eclipse), 
+0.082 to +0.09 (G400H, out-of-eclipse), and +0.02 to +0.03 (G400H, mid-eclipse). Filled squares in the G400H panel indicate the 
fluxes of the G160L order zero light computed at the same phase ranges as the G400H average spectra. Horizontal dotted lines show the 
FWHM of the G160L order zero passband. 
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Figure 2. HST lightcurves of V2051 Oph at selected passbands, for the G400H (left) and G160L (right) data. The 1421-1802 A lightcurves 
excludes the region of the CIV A1550 doublet. The curves are progressively displaced upwards by 1 mjy. Horizontal lines at mid-eclipse 
show the true zero level in each case. Vertical dotted lines mark ingress/egress phases of the white dwarf and bright spot as measured in 
section H. 
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Figure 3. Ground-based lightcurves of V2051 Oph. The white light lightcurves were artificially scaled to the same out-of-eclipse level of 
the B-band lightcurve for simplicity of display. The curves are progressively displaced upwards by 1 mjy. Horizontal lines at mid-eclipse 
show the true zero level in each case. Vertical dotted lines mark ingress/egress phases of the white dwarf as measured in section H. 
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Figure 4. Example of measuring contact phases in a lightcurve of V2051 Oph. (a) Original G160L white light lightcurvc at a time 
resolution of 5 X 10 — 5 cycles, (b) median-filtered lightcurve, shifted upwards by 0.5 mjy. (c) the derivative of lightcurve (b) after 
removing the slowly-varying disc component, multiplied by a factor 2. (d) the reconstructed white dwarf lightcurve, shifted downwards 
by 1 mjy. (e) the lightcurve without the white dwarf component, shifted downwards by 2 mjy. Dotted lines mark the contact phases of 
the white dwarf and dashed lines mark the contact phases of the bright spot. 
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Figure 5. The combined W+B ligiitcurve (runs LI to L4) and its derivative curve. Vertical dotted lines depict the median contact 
phases of the WD as listed in Tablet]. 
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Figure 6. Inferring the binary geometry from the ingress/egress phases of WD and BS. (a) Ingress-egress phases diagram. The observed 
phases of mid-ingrcss/cgress of WD are marked with crosses, those of BS with a small box. A diagonal dotted line depict the line 
joining the component stars. Theoretical gas stream trajectories for various values of q are plotted. The stream of matter passes through 
the position of BS for q = 0.19. The squashed circle represents the accretion disc whose edge passes through the position of BS, for 
q = 0.19 ±0.03. This yields a disc radius of i?(, s = 0.56 Rli- (b) the adopted geometry of the binary for q = 0.19. The observed positions 
of WD and BS are shown with the theoretical gas stream and a disc of radius Ri, a = 0.56 Rli- 
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Figure 7. Primary - Secondary mass diagram for V2051 Oph. Solid lines show the constraints obtained from the measured A w ^ and the 
inferred mass ratio of q = 0.19. Dotted lines indicate the l-cr limit on these relations. The mass function for a white dwarf radial velocity 
of K\ = 91 km/s (Watts et al. 1986) is shown as a dashed line and a mass-radius relation of R2/R0 = M2/M0 is shown as a dotted 
line. The gray cloud of points shows the confidence region and is the result of a 10 4 points Monte Carlo simulation with the values of 
A wd and q. 
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